AFRL-MN-EG-TP-2006-741 4 


A  VALIDATION  OF  FIRST-ORDER  DETONATION  SHOCK  DYNAMICS 
THEORY 


David  Lambert 

Air  Force  Research  Laboratory 
Munitions  Directorate 
AFRL/MNMW 
Eglin  AFB,  FL  32542-6810 

Sunhee  Yoo 
D.  Scott  Stewart 

University  of  Illinois  -  College  of  Engineering 
104  S.  Wright  St. 

Urbana,  Illinois  61801 

OCTOBER  2006 

CONFERENCE  PAPER 

This  paper  was  presented  at  the  2005  Winter  Simulation  Conference,  held  2-5 
December,  2005  in  Orlando,  FL.  One  of  the  authors  is  a  U.S.  Government  employee 
working  within  the  scope  of  his  position;  therefore,  the  U.S.  Government  is  joint  owner 
of  the  work.  This  paper  will  be  published  in  the  proceedings.  The  publisher  may 
assert  copyright.  If  so,  the  U.S.  Government  has  the  right  to  copy,  distribute,  and  use 
the  work  by  or  on  behalf  of  the  U.S.  Government.  Any  other  form  of  use  is  subject  to 
copyright  restrictions. 

This  paper  is  published  in  the  interest  of  the  scientific  and  technical  information  exchange. 
Publication  of  this  paper  does  not  constitute  approval  or  disapproval  of  the  ideas  or  findings. 


DISTRIBUTION  A:  Approved  for  public  release;  distribution  unlimited. 

Approval  Confirmation  #AAC/PA  06-23-05-254 


AIR  FORCE  RESEARCH  LABORATORY,  MUNITIONS  DIRECTORATE 


■  Air  Force  Materiel  Command  ■  United  States  Air  Force 


■  Eglin  Air  Force  Base 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


The  public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources, 
gathering  and  maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection 
of  information,  including  suggestions  for  reducing  the  burden,  to  Department  of  Defense,  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports 
(0704-0188),  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be 
subject  to  any  penalty  for  failing  to  comply  with  a  collection  of  information  if  it  does  not  display  a  currently  valid  OMB  control  number. 

PLEASE  DO  NOT  RETURN  YOUR  FORM  TO  THE  ABOVE  ADDRESS. 


1.  REPORT  DATE  (DD-MM-YYYY) 
XX- 10-2006 


4.  TITLE  AND  SUBTITLE 


2.  REPORT  TYPE 

CONFERENCE  PAPER 


3.  DATES  COVERED  (From  -  To) 


5a.  CONTRACT  NUMBER 


A  VALIDATION  OF  FIRST-ORDER  DETONATION  SHOCK  DYNAMICS 
THEORY 


6.  AUTHOR(S) 


5b.  GRANT  NUMBER 


5c.  PROGRAM  ELEMENT  NUMBER 


62602F 


5d.  PROJECT  NUMBER 


David  Lambert 
Sunhee  Yoo 
D.  Scott  Stewart 


5e.  TASK  NUMBER 


5f.  WORK  UNIT  NUMBER 


8.  PERFORMING  ORGANIZATION 
REPORT  NUMBER 

AFRL-MN -EG-TP-2006-74 1 4 


10.  SPONSOR/MONITOR'S  ACRONYM(S) 

AFRL-MN-EG 


11.  SPONSOR/MONITOR'S  REPORT 
NUMBER(S) 


SAME  AS  BLOCK  8 


7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory  University  of  Illinois 

Munitions  Directorate  College  of  Engineering 

AFRL/MNMW  1 04.  S .  Wright  Street 

Eglin  AFB,  FL  32542-6810  Urbana,  Illinois  61801 


9.  SPONSORING/MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES) 

Air  Force  Research  Laboratory 
Munitions  Directorate 
AFRL/MNMW 
Eglin  AFB,  FL  32542-6810 


12.  DISTRIBUTION/AVAILABILITY  STATEMENT 

DISTRIBUTION  A:  Approved  for  public  release;  distribution  unlimited. 

Approval  Confirmation  #AAC/PA  06-23-05-254 


13.  SUPPLEMENTARY  NOTES 

This  paper  was  presented  at  the  2005  Winter  Simulation  Conference,  held  2-5  December,  2005  in  Orlando,  FL.  One  of  the  authors  is  a  U.S.  Gov't 
employee  working  within  the  scope  of  his  position;  therefore,  the  U.S.  Government  is  joint  owner  of  the  work.  (*Cont.  on  abstract  below) 


14.  ABSTRACT 

*This  paper  will  be  published  in  the  proceedings.  The  publisher  may  assert  copyright.  If  so,  the  U.S.  Government  has  the  right  to  copy,  distribute, 
and  use  the  work  by  or  on  behalf  of  the  U.S.  Government.  Any  other  form  of  use  is  subject  to  copyright  restrictions. 

High  energy  explosives  are  used  in  a  variety  of  applications,  from  military  to  industrial  processes.  The  use  of  embedded,  inert  material  “wave 
shapers”  is  a  primary  method  to  customize  the  detonation  front  for  desired  ex-plosive  applications.  These  systems  create  detonation  states  that  do 
not  follow  the  simple  line  of  sight,  or  Huy-gens  model  and,  hence,  advanced  detonation  physics  with  associated  theory  are  required.  The  theory  of 
detonation  shock  dynamics  (DSD)  is  one  such  description  used  to  pro- vide  high  fidelity  modeling  of  complex  wave  structures.  A  collection  of 
experiments  using  ultra-high  speed  cameras  is  presented  as  a  means  of  obtaining  spatial  and  temporal  characteristics  of  complex  detonation  fronts 
that  validate  the  DSD  descriptions.  The  method  of  test,  operational  conditions  and  results  are  given  to  demonstrate  the  use  of  high-rate  imaging  of 
detonation  events  and  how  this  vali-dates  our  understanding  of  the  physics  and  the  capability  of  advanced  detonation  wave  tracking  models. 


15.  SUBJECT  TERMS 

Detonation  Shock  Dynamics,  Advanced  Modeling 


|  16.  SECURITY  CLASSIFICATION  OF: 

17.  LIMITATION  OF 

18.  NUMBER 

a.  REPORT 

b.  ABSTRACT 

c.  THIS  PAGE 

ABSTRACT 

OF 

PAGES 

UNCLASSIFIED 

UNCLASSIFIED 

UNCLASSIFIED 

SAR 

8 

David  Lambert 

19b.  TELEPHONE  NUMBER  (Include  area  code) 


Standard  Form  298  (Rev.  8/98) 

Prescribed  by  ANSI  Std.  Z39.18 


Proceedings  of  the  2005  Winter  Simulation  Conference 
M.E.  Kuhl,  N.M.  Steiger,  F.B.  Armstrong,  and  J.A.  Joines,  eds. 


A  VALIDATION  OF  FIRST-ORDER  DETONATION  SHOCK  DYNAMICS  THEORY 


David  Lambert 


Air  Force  Research  Laboratory,  Munitions  Directorate 
101  W.  Eglin  Blvd.,  Ste  135 
Eglin  Air  Force  Base,  FL  32542,  U.S.A. 


ABSTRACT 

High  energy  explosives  are  used  in  a  variety  of  applica¬ 
tions,  from  military  to  industrial  processes.  The  use  of 
embedded,  inert  material  “wave  shapers”  is  a  primary 
method  to  customize  the  detonation  front  for  desired  ex¬ 
plosive  applications.  These  systems  create  detonation 
states  that  do  not  follow  the  simple  line  of  sight,  or  Huy¬ 
gens  model  and,  hence,  advanced  detonation  physics  with 
associated  theory  are  required.  The  theory  of  detonation 
shock  dynamics  (DSD)  is  one  such  description  used  to  pro¬ 
vide  high  fidelity  modeling  of  complex  wave  structures.  A 
collection  of  experiments  using  ultra-high  speed  cameras  is 
presented  as  a  means  of  obtaining  spatial  and  temporal 
characteristics  of  complex  detonation  fronts  that  validate 
the  DSD  descriptions.  The  method  of  test,  operational 
conditions  and  results  are  given  to  demonstrate  the  use  of 
high-rate  imaging  of  detonation  events  and  how  this  vali¬ 
dates  our  understanding  of  the  physics  and  the  capability  of 
advanced  detonation  wave  tracking  models. 

1  BACKGROUND 

1.1  The  Need  for  Advanced  Modeling 

A  primary  driver  for  the  need  to  develop  more  efficient 
models  of  explosive  events  is  brought  about  by  the  issue  of 
length  scale.  Reaction  zones  of  detonation  fronts  are  only 
fractions  of  a  millimeter  and  using  direct  numerical  simula¬ 
tion  methods  forces  one  to  resolve  or  discretize  that  zone 
with  50  to  100  computational  “elements”  or  “cells”.  The 
shock  front  is  where  pressures  near  instantaneously  change 
from  ambient  (1-atm)  of  the  solid  phase  explosive  to  hun¬ 
dreds  of  thousands  of  atmospheres  within  the  chemical  de¬ 
composition  regime  of  the  “reaction  zone”  and  then  decay 
back  to  hundreds  of  atmospheres  of  pressure  in  the  expand¬ 
ing  gas  products  field.  State  variables  of  specific  volume 
and  temperature  also  change  rapidly  over  the  extremely 
small  reaction  zone.  The  numerical  stability  condition,  or 
Courant-Friedrich-Levy  condition,  limits  the  computa¬ 
tional  time  step  based  on  the  smallest  spatial  cell  or  ele¬ 
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ment  that  information  can  be  calculated  across.  Since  most 
munitions  are  on  the  order  of  meters,  then  these  calcula¬ 
tions  quickly  become  time  intensive  with  some  taking 
weeks  to  complete.  The  detonation  propagation  algorithms 
must  accurately  capture  the  proper  state  variables  in  order 
to  advance  the  front  and  give  accurate  time-of-arrival  and 
loading  conditions  to  the  inert  material  being  acted  upon. 
Therefore,  more  efficient  methods  are  required  to  provide 
high  fidelity  representation  of  the  physics  driven  by  small- 
scale  phenomenon  for  the  global,  larger  scale  event.  Deto¬ 
nation  Shock  Dynamics  (DSD)  is  one  such  solution  (Stew¬ 
art  and  Bdzil  1998). 

A  conventional  detonation  propagation  model  assumes 
that  the  detonation  shock  propagates  normal  to  itself  at  the 
Chapman-Jouguet  (CJ)  velocity.  This  motion  rule  is  called 
the  Huygens  construction,  or  line  of  sight.  One  can  use 
this  propagation  rule  to  propagate  a  front  through  the  ex¬ 
plosive  geometry  and  calculate  the  time  of  arrival  (TOA)  at 
the  surfaces  of  the  inert  materials.  Knowing  the  arrival 
time  of  the  fronts  from  this  simple  engineering  rule  had 
been  almost  enough  to  work  out  basic  designs.  But,  given 
the  demands  of  more  complex  systems  with  multiple  igni¬ 
tion  points,  smaller  charges,  complex  geometries  for  the 
explosive  charge,  and  the  desire  to  engineer  complex  se¬ 
quences  of  synchronous  actions,  this  simple  model  and 
motion  rule  does  not  provide  enough  accuracy. 


1.2  Review  of  Detonation  Shock  Dynamics 

The  asymptotic  theory  of  detonation  shock  dynamics  refers 
to  hydrodynamic  flow  theory  that  corrects  a  near  planar 
detonation  flow  normal  to  the  shock  to  account  for  changes 
due  to  the  shock  curvature  and  a  higher  order  theory  that 
includes  unsteady  effects.  The  theory  assumes  that  the  ra¬ 
dius  of  curvature  of  the  shock  is  large  compared  to  the 
length  of  the  reaction  zone  that  supports  the  detonation. 
The  term  "Detonation  Shock  Dynamics"  (DSD)  describes 
both  the  asymptotic  theory  and  its  engineering  application 
(Bdzil  and  Stewart  1988).  The  simplest  theory  of  DSD 
approximates  the  local  detonation  velocity  normal  to  the 
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shock,  Dn,  as  being  the  Chapman- Jouget  velocity,  DCj,  ad¬ 
justed  by  the  local  total  curvature,  /c,  relation, 


Dn=DCJ(l-vK)  (1) 

where  v  is  an  empirical  constant.  This  curvature  relation  is 
used  in  conjunction  with  an  equation  of  state  and  reaction 
rate  law  to  provide  a  new  motion  rule  that  can  be  used  to 
make  refined  engineering  prediction  of  the  detonation 
process  in  desired  applications.  The  DSD  motion  rules  are 
intended  to  replace  the  traditional  Huygen’s  construct  that 
simply  transmits  the  information  at  a  constant  DCj .  Optical 
high-speed  smear  camera  techniques  can  be  used  to  obtain 
the  Dn- k relationship  (Hull  1993  and  Lambert  2005). 

An  equation  of  state  (EOS)  is  a  description  of  the  en¬ 
ergy,  pressure,  specific  volume  and  bum  fraction  relation¬ 
ship  e(p,v,/ 1)  for  the  reacting  explosives.  The  rate  law  gov¬ 
erns  the  rate  at  which  the  reaction  progresses  from  the  solid 
phase  to  complete  gas  products  r(p,v,/ i).  The  specific 
forms  of  the  EOS  and  rate  laws  are  not  given  or  required 
here,  but  are  mentioned  for  completeness.  There  are  stan¬ 
dard  experiments,  some  using  optical  methods  (Lambert, 
2005)  for  calibrating  the  EOS  and  rate  law  models. 


2  VALIDATION  EXPERIMENTS 

A  series  of  validation  experiments  are  typically  required  to 
confirm  these  models  can  provide  the  accuracy  and  fidelity 
of  computing  the  response  in  real  explosive  applications. 
The  ability  to  accurately  track  the  wave  front  in  complex 
detonation  scenarios  is  of  high  interest  and  is  the  primary 
topic  of  interest  for  this  section.  The  Munitions  Director¬ 
ate  has  two-dimensional  (2-D)  and  three-dimensional  (3-D) 
research  codes  that  were  developed  by  the  University  of 
Illinois  specifically  for  addressing  problems  of  interest 
(Yoo  2003).  These  experiments  and  their  hardware  were 
configured  to  capture  the  physics  of  interest  within  actual 
explosive  applications  but  with  as  simplified  geometries  as 
possible  and  with  well-characterized  materials.  They  pro¬ 
vide  a  means  to  validate  the  ensemble  of  physics  models 
including  the  equation  of  state  and  rate  laws. 

2.1  2-D  Axisymmetric  “Passover”  Experiment 

One  method  of  “wave  shaping”  is  to  insert  an  inert  object 
within  the  reactive  flow  field.  This  common  method  of 
wave  shaping  causes  the  oncoming  detonation  to  track 
around  the  object  and,  hence,  allows  a  tailoring  of  its 
shape.  The  “Passover”  Experiment,  shown  in  Figure  1  was 
designed  to  provide  validation  of  a  2-D  axisymmetric  ex¬ 
plosive  design.  It  is  a  “self-light”  method  -  meaning  the 
photons  generated  within  the  detonation  reaction  zone  it¬ 
self  exposes  the  film.  The  primary  data  is  the  time-of- 


arrival  of  the  wave  exiting  the  top  surface  (explosive/water 
interface)  of  the  assembly.  The  embedded  disk  of  lead 
serves  as  an  “obstacle”  to  transform  a  simple  hemispheri¬ 
cal,  incoming  detonation  wave  to  a  reasonably  complex 
flow.  It  forces  the  simulation  tools  to  calculate  the  wave 
progression  across  a  broad  range  of  curvatures,  including 
the  transition  from  a  divergent  to  convergent  regime. 
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Figure  1 :  Cross-section  View  of  the  2-D  "Passover”  Experiment 


2.2  Experiment  Description 

The  assembly  consists  of  four  significant  components.  The 
first  is  the  main  explosive  charge  of  PBX-9501.  PBX-9501 
is  one  formulation  that  has  its  Dn-K  relationship  and  suit¬ 
able  equation  of  state/rate  law  pair  characterized.  The 
main  charge  consists  two  pieces:  a  lower  solid  cylinder  of 
69.85-mm  (2.75-inch)  diameter  by  50.8-mm  (2.00-inch) 
long  and  an  upper  cylinder  of  same  diameter  but  with 
length  20.32-mm  (0.80-inch)  and  having  a  25.4-mm  (1. 00- 
inch)  diameter  cavity  bored  to  10.16-mm  (0.40-inch)  depth 
in  one  end.  The  second  component  is  a  disk  of  pure  lead, 
machined  to  fit  the  cavity  in  the  top  explosive  piece.  The 
third  component  is  the  polycarbonate  hardware  that  physi¬ 
cally  locates  the  explosive  cylinders  together  and  holds  the 
initiation  system  precisely  at  the  centerline  of  the  base 
charge.  The  initiation  system  is  the  fourth  component  of 
the  experiment  assembly.  A  precision  initiation  coupler 
(PIC)  is  located  up  against  the  base  of  the  lower  charge.  It 
is  simply  a  steel  cylinder  with  an  "hour-glass"  internal  cav¬ 
ity  that  is  filled  with  Composition  A-5  explosive.  An  ex- 
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plosive  bridge-wire  detonator  initiates  the  PIC  that,  in  turn, 
centers  the  detonation  and  outputs  a  hemispherical  wave. 

At  time  zero,  a  single  detonator  ignites  the  bottom  sur¬ 
face.  A  hemispherical  (convex  with  positive  curvature) 
detonation  shock  front  propagates  through  the  explosive 
until  it  encounters  the  embedded,  inert  disk.  The  shock 
speed  in  the  lead  in  much  lower  than  the  detonation  veloc¬ 
ity  and,  thus,  the  wave  diffracts  around  the  disk’s  geome¬ 
try.  The  axial  symmetry  of  the  disk  and  the  encountering 
wave  transforms  the  leading  shock  into  a  torroidal  shape. 
As  the  wave  continues  around  the  leeward  side  of  the  disk 
the  “hole”  closes  up  and  interacts  in  a  symmetrical  implo¬ 
sion  type  of  event.  The  wave  structure  at  this  point  is  the 
shape  of  an  ‘m”  if  viewed  as  a  2-D  slice  through  the  cen¬ 
terline.  The  outer  regions  are  of  unperturbed,  convex  cur¬ 
vature  while  the  central  region  is  a  highly  concave  domain 
with  high  negative  curvature.  The  detonation  front  then 
reaches  the  top  surface  of  the  charge  and  is  quenched  by 
the  water. 

2.3  Diagnostics  -  The  Imacon  468  and  Cordin  132A 

Time-of-arrival  is  the  primary  data  of  interest.  The  Cordin 
132 A  smear  camera  viewed  the  top  surface  of  the  appara¬ 
tus  with  a  lOOu  wide  slit  positioned  across  the  centerline. 
The  images  in  Figure  2  show  the  smear  camera  images 
from  two  experiments.  The  camera  used  70-mm  format 
TMAX  IS03200  black  and  white  film,  writing  at 
12mm/us.  One  is  of  the  exact  configuration  of  Figure  1 
while  the  other  is  a  similar  experiment  but  done  without 
the  lead  disk  and  the  top  piece  of  explosive  (i.e.  just  the 
bottom  2-inch  length  of  explosive)  to  get  data  of  purely  a 
single  hemispherical  detonation  front.  The  figures  show 
the  wave  profile  as  it  breaks  out  on  the  explosive/water  in¬ 
terface  and  this  data  is  used  to  obtain  time-of-arrival  in¬ 
formation. 

A  DRS  Technologies  Imacon  468  equipped  with  8- 
framing  channels  was  used  in  parallel  to  the  Cordin  132A 
camera.  The  images  in  Figure  3  show  the  global  view  of 
the  event  for  the  charge  not  having  an  embedded  disk.  It  is 
simply  a  single  point  initiation  of  a  2-inch  cylinder  and  the 
images  are  of  the  detonation  front  traveling  from  the  rear 
of  the  charge  and,  in  this  vantage  point,  coming  out  of  the 
plane  of  the  page.  The  frames  are  70ns  in  exposure  with 
150ns  interframe  time.  Frames  from  1-3  show  the  glow  of 
the  hemispherical  front  within  the  PBXN-9501  charge. 
The  explosive  is  translucent  under  this  intense  light  source. 
Frame  4  is  just  after  the  hemispherical  shell  of  a  front  in¬ 
tersected  the  plane  of  the  explosive/water  interface  giving 
an  ‘O’  shape  (think  of  a  slice  through  a  hemispherical 
shell).  Frames  5-8  show  the  expanding  front  it’s  radius 
grows  and  the  truncated  hemisphere  analogy  gives  the  ap¬ 
pearance  of  a  wave  traveling  radially  outward.  The  left¬ 
most  image  of  Figure  2  is  the  simultaneous  smear  record  of 
this  event.  The  interesting  physics  aspect  of  this  image  se¬ 


quence  is  that  one  is  now  looking  at  the  instantaneous  tran¬ 
sition  of  solid  phase  explosive  to  its  gas  phase  detonation 
products.  The  imaging  time  is  below  that  at  which  the 
product  gases  have  time  to  move  outward  and  obscure  the 
view.  Pressures  on  the  outside  radius  of  the  detonation 
front  are  at  ambient,  1-atm  pressures  while  that  just  inside 
the  radius  of  the  detonation  front  are  approximately 
350,000-atm. 


No  embedded  disk  of  lead,  Embedded  disk.  FUI  apparatus 

Just  a  single  pt.  Initiation  of  2.7 5-  as  per  Figure  1 
in  diam.  By  2-inch  long  cylinder 


Figure  2:  Smear  Records  of  Passover  Test  with  No  Embedded 
Inert  (left),  with  Lead  Disk  (right) 

The  image  sequence  of  Figure  4  shows  the  wave  struc¬ 
ture  for  that  of  the  “Passover”  assembly  in  Figure  1.  The 
corresponding  smear  camera  image  is  that  of  the  right  im¬ 
age  of  Figure  2.  Frames  1-5  of  Figure  4  show  the  wave  in¬ 
side  the  charge  as  it  travels  up  and  around  the  lead  disk. 
Note  the  central  black  region  in  frames  1-4  are  because  the 
disk  is  obscuring  the  light  at  the  central  region  of  the  hemi¬ 
spherical  wave  and  the  time  of  frame  5  is  when  the  wave  is 
diffracting  around  the  disk.  Breakout  of  the  wave  on  the 
top  explosive/water  surface  occurs  between  frames  5  and 
6.  In  frame  6  one  can  see  the  torroid  shell  structure  inter¬ 
secting  the  plane  of  the  explosive  surface  and  now  there  is 
a  detonation  collapsing  radially  inward  and  one  moving 
outward.  These  opposing  detonation  fronts  continue  and 
frame  8  is  opportunistic  because  it  captures  the  collision  of 
the  inward  running  wave  as  it  is  just  collapses  on  center- 
line.  These  images  were  captured  using  60ns  exposure  du¬ 
ration  with  150ns  constant  interframe  time. 
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Figure  3:  Imacon  468  Image  Sequence  of  Single  Point  Detona¬ 
tion  of  a  Right  Circular  Cylinder  of  PBX-9501  (No  Inert  Disk) 


Figure  4:  Imacon  468  Sequence  of  the  “Passover”  Experiment 
using  Hardware  in  Figure  1. 
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2.4  Results  and  Comparison  With  DSD  Theory 

Four  experiments  of  the  Passover  configuration  were  con¬ 
ducted.  Their  smear  records  (all  were  similar  to  that  of  the 
right  image  in  Figure  2)  were  digitized  to  give  time-of- 
arrival  across  the  diameter  of  the  charge.  This  data  is 
shown  in  Figure  5  along  with  that  from  a  2-D  simulation 
code,  called  “2-D  Wavetracker”,  that  has  the  DSD-based 
motion  rule.  Also,  superposed  is  the  Passover  experiment 
simulated  using  common  Huy  gen’s  construct  and  that  of  a 
direct  numerical  simulation  (DNS).  The  DNS  is  consid¬ 
ered  an  exact  solution  by  solving  the  full  set  of  Euler  equa¬ 
tions  with  high-order  solvers.  Recall  though,  that  the  DNS 
represents  the  computationally  intensive  method  that  the 
DSD-model  seeks  to  relieve.  All  data  is  given  with  the 
first  breakout  time  adjusted  to  give  the  collection  the  best 
overall  fit. 

The  DSD  and  DNS  records  are  nearly  the  same  with 
an  absolute  time  difference  of  less  than  0.05us,  except  in 
the  center  region..  The  Huygen’s  construction  assumed 
constant  normal  shock  speed  of  DCJ  =8.86  mm/us.  The 
time-of-arrival  at  the  top  surface  of  the  Huygens  is  about 
0.200us  shorter  (overall)  than  that  of  either  the  DSD  or 
DNS  results.  The  DSD  result  has  greater  discrepancy  with 
the  DNS  in  the  central,  colliding  wave  region  that  may  in¬ 
dicate  inappropriateness  of  the  first  order  approximation 
and  its  assumption  of  small  reaction  zone  thickness-to- 
local  curvature. 

k  -  —  «  1  where  lr  =  reaction  zone  thickness  (2) 


A  2-D  axisymetric  view  of  the  isochronal  plot  for  the 
DSD-model  is  given  in  Figure  6.  One  can  see  how  the 
wave  front  is  transformed  from  the  hemispherical  shell  to 
the  torroidal  shape  as  it  encounters  and  sweeps  the  lead 
disk  to  the  interacting  structure.  Note  the  variation  in  cur¬ 
vature  of  the  resulting  wave  structure  as  it  exits  the  top  of 
the  assembly. 
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Figure  5:  Passover  Experiments  (“PASS-1”,  “PASS-2”,  “PASS- 
3”,  and  “PASS-5”)  Compared  with  DSD-model  Motion  Rule, 
Huygens  Wave,  and  Direct  Numerical  Simulation  (DNS). 


Figure  6:  Isochronal  Contours  for  the  DSD-based  Motion  Rule 
Model  of  the  Passover  Experiment. 


3  CONCLUSIONS 

The  Passover  experiment  generates  a  large  variation  in  the 
normal  detonation  velocity  and  produces  a  reproducible 
and  tractable  means  by  which  to  validate  detonation  front 
tracking  models.  The  selection  of  time-of-arrival  of  the 
shock  at  a  known  location  is  but  one  of  many  parameters 
of  possible  comparison.  For  many  explosive-metal  appli¬ 
cations,  it  is  a  critical  parameter  and  should  be  considered 
a  fundamental  one  that  simulations  should  match.  The 
Passover  configuration  represents  such  applications  but  in 
a  well-controlled,  well-characterized  experimental  test  bed. 

The  DSD-WaveTracker  tool  with  its  associated  equa¬ 
tion  of  state  and  rate  law  shows  promise  as  a  computation¬ 
ally  efficient,  but  accurate,  model  to  use  in  further  designs 
and  applications.  The  absolute  difference  of  0.04us  in 
TOA  between  the  DSD  and  DNS  results  is  only  a  0.4%  dif¬ 
ference  for  the  total  time  of  wave  propagation  along  the 
length  of  the  charge.  The  Huygen’s  wave  construct  (recall 
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that  this  is  the  predominately  used  rule)  was  shown  to  be 
an  order  of  magnitude  larger  in  time  discrepancy  of  the 
waves  arrival. 

Excellent  agreement  between  experiment,  theory,  and 
direct  numerical  simulation  was  achieved.  The  level  of 
agreement  is  encouraging  because  it  indicates  that  one  will 
be  able  to  use  these  models  for  real  explosives  and  more 
accurately  predict  the  response  of  complex  systems. 

REFERENCES 

Bdzil,  J.B  and  D.S.  Stewart.  1989.  Modeling  of  two- 
dimensional  detonation  with  detonation  shock  dynam¬ 
ics.  Physics  of  Fluids  A,  1:  1261-1267. 

Hull,  L.  1993.  Mach  reflections  of  spherical  detonation 
waves.  In  Proceedings  of  the  10th  International  Sym¬ 
posium  on  Detonation. 

Lambert,  D.E.  2005.  Experimental  validation  of  advanced 
explosive/metal  interactions.  In  Proceedings  of  the 
Conference  on  High-Speed  Photography  and  Photon¬ 
ics,  SPIE  International  Symposium  on  Optics  and 
Photonics. 

Stewart,  D.S.  and  J.B  .  Bdzil.  1988.  A  lecture  on  detona¬ 
tion  shock  dynamics.  Mathematical  Modeling  in 
Combustion  Science,  Lecture  Notes  in  Physics  249: 
Springer- Verlag,  17-30. 

Yoo,  S.  2003.  Wave  tracking  in  complex  geometries.  Doc¬ 
toral  dissertation,  Theoretical  and  Applied  Mechanics, 
University  of  Illinois. 


currently  working  as  a  Research  Engineer  in  D.  S.  Stew¬ 
art's  research  group  in  the  Department  of  Mechanical  and 
Industrial  Engineering  at  the  University  of  Illinois.. 

D.  SCOTT  STEWART  obtained  a  BS.  E.  at  the  Univer¬ 
sity  of  New  York  at  Buffalo,  Engineering  Science  in  1976 
and  his  Ph.  D.  in  Theoretical  and  Applied  Mechanics,  Cor¬ 
nell  University  in  1981.  He  joined  Theoretical  and  Ap¬ 
plied,  University  of  Illinois  in  1981  and  is  a  full  professor. 
His  work  has  been  generously  supported  by  Los  Alamos 
National  Laboratory  and  the  U.  S.  Air  Force  in  the  areas  of 
reactive  flow,  shock  physics,  modeling,  advanced  simula¬ 
tion  and  miniaturization  of  explosive  technology.  He  is  a 
fellow  of  both  the  American  Physical  Society,  the  Institute 
of  Physics  and  an  Associate  Fellow  of  AIAA.  He  is  cur¬ 
rently  a  Professor  of  Mechanical  and  Industrial  Engineer¬ 
ing  at  the  University  of  Illinois. 


AUTHOR  BIOGRAPHIES 

DAVID  LAMBERT  is  the  Senior  Scientist  in  the  Damage 
Mechanisms  Branch  of  the  Munitions  Directorate,  Air 
Force  Research  Laboratory.  His  research  interests  are  in 
the  area  of  detonation  physics  and  high-rate  material  be¬ 
havior.  He  received  his  BS  in  Mechanical  Engineering 
from  Florida  State  University,  MS  in  Engineering  Mechan¬ 
ics  and  PhD  from  the  University  of  Florida.  He  is  a  mem¬ 
ber  of  APS,  AIAA,  and  NDIA-National  Defense  Industrial 
Association. 

SUNHEE  YOO  obtained  his  B.  A.  in  Mathematics  from 
Kyung-Pook  National  University,  Taegu,  Korea  in  1987, 
M.  S.  in  Applied  Mathematics  from  Korea  Advanced  In¬ 
stitute  of  Science  and  Technology,  Seoul,  Korea  in  1989. 
At  the  University  of  Illinois  he  received  a  M.  S.  in  Mathe¬ 
matics  in  1998  and  a  Ph.  D.  from  Theoretical  and  Applied 
Mechanics  in  2003.  While  in  Korea  he  worked  in  industry 
for  Hyundai  Electronics  as  a  software  engineer,  and  later  in 
Champaign,  Illinois  for  the  U.  S.  Army  Construction  Re¬ 
search  Laboratory  from  1994-1988  while  studying  mathe¬ 
matics.  His  current  interests  are  in  computational  mathe¬ 
matics,  simulation  technologies  and  shock  physics,  and  is 


